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We show that, at low temperatures, macroscopic inhomogeneities of the elec- 
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■ tron density in the interior of a finite sample cause a reduction in the measured 
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conductivity peak heights cr^ ax compared to the universal values previously 
predicted for infinite homogeneous samples. This effect is expected to occur 

in ; 

ON ■ for the conductivity peaks measured in standard experimental geometries such 



as the Hall bar and the Corbino disc. At the lowest temperatures, the de- 
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*^ ■ crease in cr x c x 3X (T) is found to saturate at values proportional to the difference 

a : 

between the adjacent plateaus in a xy , with a prefactor which depends on the 



particular realization of disorder in the sample. We argue that this provides 
a possible explanation of the "non-universal scaling" of cr^ ax observed in a 
number of experiments. We also predict an enhancement of the "non-local" 
resistance due to the macroscopic inhomogeneities. We argue that, in the Hall 
bar with a sharp edge, the enhanced "non-local" resistance and the size cor- 
rections to the "local" resistance R xx are directly related. Using this relation, 
we suggest a method by which the finite-size corrections may be eliminated 

from R xx and R xy in this case. 
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I. INTRODUCTION 



The fascinating property of the quantum Hall effect (QHE) that initially attracted such 
great attention to the phenomenon is the precise quantization of the Hall conductivity o xy at 
certain values of the magnetic field. Most theoretical research has focused on the properties of 
the electron system inside these quantized plateaus in a xy . The plateaus have been associated 
with the incompressibility of the 2D electron gas, arising either from Landau quantization, 
at integer filling factors, or from electron-electron interactions, at fractional filling factors. 
The transition regions, where a xy crosses over between quantized values and the longitudinal 
conductivity a xx experiences maxima, have received less attention. The main factor which 
inhibits progress in this direction is the lack of reproducible experimental results on the 
inter-plateau regions, despite the impressive stock of data on the QHE which has been 
accumulated over the last decade. In addition to the fact that the general behavior of the 
QHE depends on the electron density, temperature and disorder, samples cut from the same 
substrate and measured at the same temperature often reveal different dependences a xy (B) 
and a xx (B). This annoying data dispersion is particularly apparent at low temperatures. 
Certain success in obtaining reproducible data has been achieved only for the critical 
behavior of the width of inter-plateau regions at low temperatures. ]2|] However, as far as 
the heights or shapes of the peaks in a xx are concerned, the general impression is that too 
many factors are involved to allow any systematic conclusions. 

On the other hand, there do exist a number of theoretical works which argue that certain 
universal behavior of the conductivity tensor must exist at low temperatures in the regions 
between well-pronounced pairs of plateaus (critical regime). Kucera and Streda H con- 
sidered semi-classical single-electron transport in a partially filled Landau level for a simple 
model of a periodic long-range potential. They found that the maximum value of a xx reached 
at half-integer filling factors does not depend on either the magnitude of the potential or 
the Landau level number N and is equal to e 2 /2h. This result was later mapped onto the 
fractional regime by using two related approximations of the correlated electron state: the 
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dirty boson j|] and the composite fermion approach. || A different sort of argument for both 
integer and fractional regime was presented in Refs. || and J7|. Based on the rather general 
assumption that, at low temperatures, the electron system in a critical transition region can 
be represented by a random mixture of two quantum liquids with different quantized local 
Hall conductivities a 2 and ax, it was shown that a xx and connected by a universal 

relation. The peak height cx^ ax was found to be equal to one half of the difference between 
the Hall conductivities of the adjacent plateaus |cr 2 — <T\\/2. For the integer peaks, this 
result yields the value e 2 /2h obtained in Ref. ||. For the fractional regime, it matches the 
results of Refs. 0] and ||, after the latter are somewhat corrected to allow for the fact that 
the maxima of a xx (B) of the integer peaks do not map exactly onto the maxima of a xx (B) 
for the principal series. || We can also refer the reader to quantum Monte-Carlo studies in 
Ref. || where the same value, 0.5e 2 /h, for the integer peaks was obtained in a simulation of 
single-electron scattering off short-range impurities. Thus, while different theoretical models 
agree on an expected universality of the conductivity peak heights and even on their values, 
experiment offers no evidence to support this prediction. 

To make the situation even more confusing, a puzzling feature was observed in a number 
of experiments performed at very low temperatures (15-40 mK): the relative heights of 
most of the conductivity peaks obtained in the fractional regime were, indeed, found to 
scale approximately as \a 2 — o~i|, but with an absolute factor which differed from 1/2 and, 
moreover, varied from sample to sample. [HJ Such a "universality within one sample", while 



there is none between different samples, is hard to understand. This feature was especially 



well seen in recent experiments performed in the Corbino geometry [TT| where scaling of the 
peak heights was observed simultaneously for the fractional and integer regimes, i.e., in quite 
different ranges of magnetic field. As a function of temperature, the height of each peak was 
found to pass through a maximum value. These maximum values differed for different peaks, 
being scattered below, though not very far from, the universal values predicted by theory. 
On the low temperature side of these maxima, therefore, with decreasing temperature the 
peak heights fell further below the universal values, as has also been earlier observed in 



low-mobility samples. Jl| At the lowest temperatures (14 mK), however, different peaks in 
a given sample converged to the corresponding theoretical values multiplied by the same 
sample-dependent factor. 

In this work, we suggest an explanation of the phenomenon of the "non-universal" scaling 
of the peak heights. Our explanation assumes that in the interior of the sample there exist 
random inhomogeneities of the electron density with a very large correlation length. While 
the origin of these inhomogeneities is unclear, we note that the existence of such fluctuations 



was also a necessary assumption in recent work [12| which proposed an explanation for 
another experimental puzzle - the Resistivity Law, which is observed in some samples at 
higher temperatures. We can only speculate as to whether such fluctuations might result 
from imperfections in the doping process or from incomplete equilibration of the electron 
density on sample cool-down, or from both. The fact that the experimental traces in the 
inter-plateau regions often change after the sample is reheated and cooled down again, 
indicates that that the latter mechanism may be important. By whatever means they might 
arise, we argue that even small fluctuations of the electron density become crucial at low 
enough temperatures, where the dependence of the conductivity tensor on the local value of 
the filling factor is almost singular. In particular, we show that, due to the finite size of the 
sample, the sample edge gives rise to a positive contribution to the Corbino resistance which 
is proportional to 1/(02 — o^). This contribution increases with decreasing temperature 
due to the growing correlation radius of the random clusters which are responsible for the 
current transfer. At sufficiently low temperatures, when the correlation radius exceeds the 
sample size, the edge contribution dominates the sample conductance. The heights of the 
observed peaks in a xx in a given sample are found to freeze at values which differ from the 
corresponding "universal" values in an infinite sample by a random geometric factor which 
is the same for all peaks. 

By making use of a current-voltage duality that exists in two-dimensional conductivity 
problems, we show how our results may also be applied to the Hall bar geometry. We find 
that macroscopic inhomogeneities have the same effect on the measured peaks in p xx in this 



geometry as on the peaks in a xx for the Corbino disc. This is in agreement with the fact that 



sample-dependent scaling of the peak heights has been observed in both geometries ||10| , |lT 

It has previously been realized that density inhomogeneities may distort conductivity 
measurements in the Hall bar geometry. In particular, for certain gated or mesa-etched 
systems, the boundary to a vacuum is believed to involve a rather gradual decrease in 
electron density. This creates a strip along the edge with a quantized Hall conductivity and 
very low scattering, while the bulk of sample can be in a transition region with noticeable 
dissipation. Such a strip can trap a significant portion of the current, making the current 
distribution in the sample inhomogeneous and affecting the measured resistance R xx . In 
the language of the edge-transport theory, this can be re-formulated as poor equilibration 
between the edge states and the bulk. While this model has been successful in accounting 
for certain non-local resistance measurements, the presence of a smooth edge cannot explain 
the observed non-universal scaling of the conductivity peak heights. Therefore, since the 
random bulk inhomogeneities are crucial for this scaling, and since the effects of a smooth 



edge have been discussed before, [1^-13 in this work we will focus only on the consequences 



of macroscopic bulk inhomogeneities. The results we present for the Hall bar therefore 
apply to samples with sharp edges (i.e., with the edge width less than magnetic length). For 
certain samples, some combination of the two models may be appropriate. 

As well as giving rise to "non-universal scaling" of the peak heights at the lowest temper- 
atures, we show that macroscopic inhomogeneities in the bulk also lead to an enhancement of 
the non-local resistance. In fact, by developing a general "boundary-strip" approach which 
describes finite-size effects in a macroscopically inhomogeneous sample, we show that, for a 
sample with a sharp edge, both of these effects are directly related. Using this relation, we 
propose an experimental method by which the edge contribution in the observed resistivity 
may be separated from the bulk contribution, at least over the range of temperature for 
which the edge contribution is not too big (at intermediate T). 

The paper is organized as follows. In Sec. II, we discuss the role of macroscopic inho- 
mogeneities in the transport properties of an infinite QHE system in different regimes of 



temperature. Finite-size effects on the low-temperature two-terminal resistance of an in- 
homogeneous Corbino disc are considered in Sec. Ill and in the Appendix. The boundary 
impedance matrix formalism, which provides a description of small edge contributions to 
the measured resistances in any sample geometry (Corbino disc, Hall bar etc.), is developed 
in Sec. IV. The explicit form of the impedance matrix is evaluated for a macroscopically 
inhomogeneous sample with a sharp edge. The enhanced "non-local" resistance arising from 
the macroscopic inhomogeneities is studied in Sec. V, where it is also shown that this is 
directly related to the edge corrections in the measured "local" and Hall resistances. Sec. 
VI concludes the paper. 

II. INFINITE SAMPLE 

We consider an infinite sample in the presence of macroscopic inhomogeneities of the 
electron density. The local value of the filling factor can be written as u(r) = V + 5v(r), 
where V is the average filling factor, and bvir) is a small fluctuating component with a 
magnitude Svq <C V. We assume that some impurity scattering occurs on scales much 
smaller than R c , such that the local conductivity tensor exhibits the quantum Hall effect. 
Specifically, we will assume that at a critical value z/ c , the Hall conductivity undergoes a 
sharp crossover between two quantized values, a± = e 2 i>i/h, and a 2 = e 2 V2/h. Here v\ 
and z/2 are a pair of adjacent values of v at which the electron system is incompressible, 
v\ < v c < V2- The diagonal conductivity (J X x[y) has a sharp peak in the crossover region, 
being very small everywhere else. The width of the crossover region 8vt vanishes as T — > 0. 
We can consider this picture as simply following from numerous experimental data on the 
low-temperature quantum Hall effect in large samples. 

The critical behaviour of the conductivity tensor has been much studied. The width of 
the crossover region vanishes as T — > according to 

with an exponent k which has been argued to take the universal value 3/7. [16| In the 
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zero temperature limit, the maximum value of the dissipative conductivity is thought to 
approach e 2 /2h in the spin-split integer quantum Hall effect, || and (o 2 — 0"i)/2 for the 
fractional regime. |6|,[7| As explained in the introduction, this prediction was deduced from 
the hypothesis that the correlated electron state in the vicinity of a QHE transition represents 
a random mixture of two incompressible liquids with localized quasiparticles on top. The 
correlation radius R™ 1C of this two-phase system was assumed to be larger than the magnetic 
length (hc/eB) 1 / 2 . To avoid confusion, we emphasize that the correlation radius of the 
macroscopic density fluctuations R c that we consider in the present work is assumed to 
be much larger than Rf lc which is thus considered as a microscopic length. The "local" 
conductivity tensor a{r) introduced above is defined at scales larger than i?™ lc (and, if 
effects of quantum interference are important, larger than the phase-breaking length) but 
smaller than R c . For our purposes, the peak value of the local dissipative conductivity will 
not be important; we shall only assume that it is either of the order of or less than the 
difference in the quantized values of the Hall conductivity, 02 — 01. We believe that this 
assumption agrees with existing data on the quantum Hall effect in a strong magnetic field. 

Let the magnetic field be tuned so that the average filling factor V is close to u c . At 
sufficiently low temperatures, the crossover width 5vt becomes much less than the fluctuation 
magnitude 5uq which is small but temperature independent. Hence, as temperature goes 
down, the conductivity distribution becomes strongly inhomogeneous. In most of the sample, 
as illustrated in Fig. l(a,b), the Hall conductivity a xy is quantized at either 01 or 02, and 
the diagonal conductivity a xx is very small. Only a narrow intermediate region within the 
interval \u — u c \ ~ Svt, has a noticeable a xx (this region is indicated by the grey color in Fig. 
1(c)). As long as V stays within this interval, the "grey" region forms an infinite percolation 
cluster. The latter consists of strips of a small width wt ~ R c &vtI& v q which join near saddle 
points (critical vertices). The characteristic size of a cluster cell, i.e., the distance between 
two vertices follows from classical percolation theory [T7J 



£ T ~ R c {5^/5u T ) A '\ (1) 
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We note that the true geometry of the percolating cluster includes also loops and by-passes 
on links between vertices not shown in Fig. 1(c). 

The evaluation of the effective conductivity tensor o* requires an understanding of how 
the current density is distributed in such an inhomogeneous system. The current distribution 
is known to be quite different in the two limiting cases which are reached as the relative 
values of the dimensionless parameters wt/R c and c^ ax /(cr 2 — 01) are varied. 

When a xx (r) is very small (how small will be determined below), the continuity con- 
ditions force the current to flow almost exactly along lines of constant a xy . Jl8|,|l5j For an 
infinite random system, this means that the current is concentrated near the percolation 
threshold of the function a xy (r), within a sparse percolation network inside the grey area 
shown in Fig. 1(c). The characteristic parameters of this new network - the width w and 
the length / of an elementary link of the cluster - are determined by the value of a xx at the 
percolation level in v. They can be estimated from the self-consistent condition that the 
current is able to cross the lines of constant a xy in order to pass from one critical saddle 
point of the network to the next one, which has a slightly different value of a xy . |T7|'P3 We 
give here the final expressions for both parameters 

I v„. V' a /RA 10/13 

V cr xx w T J 

The net diagonal conductivity of the system a* x is determined by the geometry of an ele- 
mentary link, as given by a* x ~ u xx l/w which yields, for the maximum value of a xx 

I r \ 10/13 

a:r x ~(<r) 3/13 (^-^i) 10/13 (g) • (3) 

Note that the effective conductivity decreases with increasing temperature, more or less as 
1/5ut, unless the temperature dependence er^Jf^T) is very sharp. 

The net Hall conductivity a* y does not depend on the geometry of the network and 
simply coincides, apart from small corrections, with the percolation threshold value <J xy (v). 
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Hence, a* , crosses over from o\ and a? within the same interval of v as does the local Hall 
conductivity, \V — u c \ ~ Svt- 

This picture is correct provided a xx is sufficiently small and the grey regions are not too 
narrow, such that the current can stay within these regions, w <C Wt- The last condition 
can be written as T ^> T s \, with T s % given by the equation 



The regime exists only if a xx is much less than 02 — <J\ at temperatures at which the system 
can already be considered to be macroscopically homogeneous, 8ut(T) ~ 5u Q . 

As the temperature is reduced below T sl , the current spills out into the quantum Hall 
regions. In the limit T <C T s i, which will be mostly considered in the rest of the paper, 
the currents flow predominately in the quantized Hall regions. The details of the grey areas 
become unimportant (except near the saddle points), and they may be replaced by sharp 
boundaries. This represents a particular example from a class of "black-and-white" systems 
which was studied in Refs. |!9|j6|,|7j . Due to the continuity conditions and the absence of 
scattering in the bulk, the currents in the "black" and "white" phases cannot cross the 
phase boundary. Instead, the currents have to focus at vertices to pass between the two 
corners of the same color. The net conductivity tensor depends on the ratio of the average 
currents flowing in the two phases which is controlled by V. At v c —V ^> 5ut, all saddle points 
have values of a xy close to <j\ which ensures good percolation in the "black" color, so that 
white current is relatively small and a* y ~ o\. At V — v c ^> 5ut, we get good percolation in 
"white", and a* y ~ 02- A crossover takes place within the interval 5vt when both currents 
are comparable. Since the local diagonal conductivity is small in both quantized regions, 
the net conductivity <r* x is also small when either of these phases percolates freely, and 
experiences maximum in the crossover region. In the limit of zero dissipation in each of the 
quantized Hall regions, the dependence of a* x vs. a xy of the isotropic two-phase system is 
known to be a universal function - a semicircle [0] 




(4) 




(5) 



9 



In particular, the maximum value of the average conductivity is 

eC" = (*a " *i)/2, (6) 

independent of the details of the "grey" areas. 

To conclude, the long-range inhomogeneities in the infinite sample do not essentially 
alter the conductivity peak width. They cause the peak height in a* x to saturate at low 
temperatures at the universal value (||). If the microscopic conductivity tensor components 
CxxyVxy, as argued in Refs. |||7|], also satisfy the "semicircle" relation then the presence 
of the macroscopic inhomogeneities has no effect in the low-temperature limit. In a finite 
sample, however, as we now proceed to show, a* x at T = deviates drastically from the 
universal value due to the inhomogeneities. 



III. FINITE-SIZE EFFECTS IN THE CORBINO GEOMETRY 

To study finite-size effects, we have to resort to a realistic experimental set-up which 
means that we have to specify the geometry of the sample and the attachment of the contacts. 
Experiments are performed on the Hall bar, in the van der Pauw method, and on the Corbino 
disk. Our initial choice will be the Corbino disc on which recent experiments in Ref. ITT 



have been performed. In Sec. Ill, we show how to transfer our results to the Hall bar using 
the current-voltage duality. We will not discuss the van der Pauw geometry explicitly in 
this work. 



A. Contact resistance 

A Corbino disc cut from an inhomogeneous sample is shown schematically in Fig. 2(a). 
The two-terminal resistance between the metal probes attached at the inner and outer 
circular edges is measured. We consider an ideal contact without any tunnel barriers or 
dielectric layers between metal and sample, the metal boundary having a constant potential. 
We assume, first, that the correlation radius of the "grey" cluster £j is much less than all 
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sample dimensions, ri,T2 and W = r2 — r±. In the limit £,t/W 0, the sample can be 
considered to be homogeneous with the conductivity tensor a*. The two-terminal resistance 
is then given by 

Ro = ^, A = ±\n r A (7) 

a* xx 2tt n 

where is the geometric aspect ratio, which for W <C rx,r2 is close to Aq ~ WjlitT. 
A finite value of ^t/W , as we now demonstrate, leads to an increase in the two-terminal 
resistance above this value (0) due to an effective contact contribution. 

The simplest way to estimate this correction is by monitoring the Joule heat dissipated 
in the sample. At low temperatures, T <C T sl , as already discussed, the currents flow in 
the (black or white) quantized Hall regions. Since scattering in the bulk is negligible, the 
currents cannot cross the phase boundaries and , in order to pass through the sample, must 
focus at vertices formed by adjacent corners of different phases. All the energy dissipation 
in the system occurs in "hot spots" at the vertices. The vertices are of two kinds: 4- vertices 
in the sample interior, formed by alternating "black" and "white" corners, and 3-vertices at 
the boundary, "black-white-metal". The above resistance (|7|) corresponds to dissipation at 
internal vertices, as given by 

Qint = I 2 Rq- (8) 

Since the metal contact boundary is at a constant potential, the current lines have to focus 
at 3-vertices in order to enter the (e.g., inner) contact. This causes additional dissipation at 
the edge given by 

Qedg = ^ I n Rg, 

a 

where is the current passing to the metal at the 3-vertex a, and R a is the effective 
resistance of the vertex. 

It turns out that focusing of the current lines occurs only at those of the 3-vertices 
at which the Hall conductivity a xy increases from a\ to 02 to the right when looking into 
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the sample from the metal (we call these "active" vertices). For the remaining half of the 
vertices, I a = 0. As shown in the Appendix, the resistance of an active 3-vertex is always 
given by 

R a = R 3 = — — -, (9) 

2{a 2 - ax) 

independent of the microscopic details of the vertex core. In particular, the result remains 
valid if the phase boundary branches when approaching the metal forming a fork-like struc- 
ture, as illustrated in Fig. 2(b) in magnified view. Branching, though not shown in Fig. 
2(a), does occur in a random percolation cluster on scales smaller than Only the fork 
vertices with white as the rightmost color are active. 

For a uniform pattern of white/black regions, the current entering the contact will be 
shared approximately equally between the titx/^t vertices at the edge so that 

Q% ~ ^3^. (io) 

The total two-terminal resistance _R t ot can be found from the total Joule heat, 

Qtot = I 2 Rtot = Qedg Qedg Qint- 

Using Eqs. (0H0), we find 

a* x \n r 2 / 27i{a 2 - (Jx) 

The second term in this expression is, of course, an estimate. In the next Section, we re- 
derive Eq. flTTl ) rigorously for a simple model of a periodic chessboard two-phase distribution 
with £ T replaced by the square size. 

Thus, at the maximum of (T^. x {v), the relative correction to the peak resistance arising 
from the edges is of the order of £,t/W . (We assumed here that r% is not much less than r 2 , 
and used the fact that, at the lowest temperatures, cr*™ ax ~ o 2 — ax). The experimentally 
measured value of the diagonal conductivity <r^ p defined by |H| 

a Tx = A / i?tot 
12 



is now lower than the "bulk" conductivity a* x by the same relative amount. This negative 
correction increases as temperature is lowered, since £t becomes larger [Eq. ([I])]. We believe 
that the decrease in cx^ p that we predict from these considerations provides a possible 



explanation for the observations of Rokhinson et al. [[Tl|] . 

Note that the correction to cr^f may be different for different conductivity peaks since 
the width 5vt and, hence, the correlation length £r, may vary from peak to peak. This 
accounts for the dispersion of traces 0"^ ax (T) observed for different peaks in the integer 
regime, see data for two samples by Rokhinson et al. [11] which we reproduced here in 
Fig. 3. At the lowest temperatures, however, as one can see from the figure, all the traces 
converge and tend to collapse onto one curve. The corresponding low-temperature value 
can be written as cr^ ax = ke 2 /2h, where the coefficient k is almost the same for different 
peaks, but is sample-dependent |L1| (compare Figs. 3(a) and 3(b)). As found in the quoted 
work, the well-pronounced peaks in the fractional regime have the height A; (02 — o"i)/2, 
with approximately the same k as do the peaks obtained on the same sample in the integer 
regime. This agrees with earlier data on the fractional quantum Hall effect mostly obtained 



in the Hall bar geometry [TO] which show that, at the lowest temperatures, most of the peak 
heights are proportional to the difference in adjacent plateaus in a xy , with a prefactor that 
fluctuates from sample to sample. To understand the origin of the curious "universality 
within one sample" let us consider the lowest temperatures. 



B. Low-temperature limit 

Since the characteristic correlation radius £t grows as the temperature is reduced, it 
will eventually become of the same order as the distance between the contacts W at some 
temperature T ~ T s2 . At this point, the edge contribution to the resistance in Eq. ( [TT]) 
is as large as the bulk contribution. When temperature is decreased further, T <C T S 2, a 
strong inequality £t 3> W is met for any peak, so that no critical 4-vertices (shown in Fig. 
2(a) for higher temperatures) fall within the sample. The current is transported directly 
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from contact to contact by one or more pairs of "white" and "black" clusters connecting the 
contacts. "Black" and "white" clusters of the size W can exist simultaneously while V is in 
the interval \V — u c \ ^5 5u w , where 5v w is determined from the equation 

W-Rtp-Y 3 . (12) 

Outside of this interval, percolation between contacts can exist only in one phase, either in 
the "white" or in the "black" . The resistance in this case is equal to that of a homogeneous 
Corbino sample with no scattering inside, that is, R to t = oo, and cr^ p = 0. Hence, 5u w in 
Eq. fll2|) represents the observed peak width. 

Notice that Su w does not depend on temperature and is much larger than the value for 
the infinite sample Sut- Since no 4-vertices are now found in the sample, the finite value 
of Svt is no longer relevant. All temperature dependence has now disappeared from the 
problem and all the peaks (J^P(v) must become identical. We suggest that this explains 
why the experimental traces u^iT) in integer regime converge at low temperatures |p~Q| , [TTj] . 
Recall that it was the differing values of 5vt that caused the contact contributions to the 
total resistance to vary for different peaks at T ^> T s2 . 

This saturation of the peak width which we predict at low temperatures is similar to the 
well-known saturation effect which is expected to occur in small samples when the coherence 
length £ c becomes comparable to the sample size. |21] The observable difference between the 



two mechanisms is in the size dependence, 5v w oc W K : we predict the classical index 
k = 3/4 instead of the 3/7 which is thought to be appropriate for the quantum problem. 

We now discuss the height and possible shape of the peak. As we shall show below, in 
the zero-temperature limit, the two-terminal conductance 1 / -Rtot can take only quantized 
values 

l/R tot = M{a 2 -(J 1 ), (13) 

where M is an integer (including zero) which depends on the specific realization of the 
disorder. As the average filling fraction is varied, the value of the integer M may change, 
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and the peak in cr^P(v) can display an unusual step-like dependence on filling fraction, as 
illustrated schematically in Fig. 4c and 4d. Although the shape of these peaks (number and 
position of the steps) depends on the specific realization of disorder, this shape is the same, 
with the peak height expressed in units of 02 — 01, for all peaks in a given sample. 

Let us fix the average filling factor V somewhere within the peak width 5u w , e.g., at the 
point V = v c . We have to distinguish two major cases depending on the random configu- 
ration of a sample. Either both "black" and "white" percolate in the radial direction, or 
both "black" and "white" percolate in the azimuthal direction, as illustrated in Fig. 4(a) 
and 4(b). Together with the case in which only one phase percolates in both directions, 
which we determined above as being outside of the peak width, this exhausts all topological 
possibilities. In the case of azimuthal percolation, Fig. 4(b), the current can barely pass 
between contacts since it is not allowed to cross the phase boundary. Hence, peaks are 
missing, R tot = 00. More precisely, a^P is as small as the average of a xx in the quantized 
regions, and therefore rapidly vanishes as T — > 0. Such behavior is observed in experiment 
25], but it is usually attributed to a "bad sample" or "bad contacts" and, consequently, 



does not reach publication. In some sense, it is correct to say this in our model as well; 
however, both "bad" and "good" samples belong to the same statistical ensemble with a 
quite small amplitude of inhomogeneities. In the case shown in Fig. 4(a), when the contacts 
are connected by one "black" and by one "white" region the resistance is finite and equal to 
the doubled resistance of the 3- vertex i? tot = 2R 3 = l/(a 2 — see Eq. @ for the 3- vertex 
resistance. 

We can easily obtain the shape of the peak in a^P for the configuration in Fig. 4(a). Let 
v\ and V2 be the filling factors at the two saddle points which control the current transfer, 
v-i — V\ ~ 8v w . When V is shifted down and crosses the level of the lower saddle-point z/ 1; the 
corresponding "black" bridge becomes "white" , and percolation in "black" quits, so that the 
resistance i? to t becomes infinite. This transition occurs abruptly in V (more precisely, within 
a small interval ~ 8vt)- Similarly, for v > v 2 the "white" region percolates, and the system 
is on a quantized Hall plateau corresponding to z/ 2 . The resulting "peak" in <r™ p represents 

15 



a box of a width V2 — V\, Fig. 4(c). The latter is of the order of 8v w and fluctuates from 
sample to sample by 100%. 

There may, of course, be more than one pair of "black" and "white" regions connecting 
the contacts (and, hence, more saddle points which switch). Since all of the 3- vertices at the 
same edge are connected in parallel, the inverse resistance is given by Eq. (|T3"D, where M is 
the number of pairs of connecting regions. In the degenerate case shown in Fig. 4(b), M = 0. 
Thus the total conductance of the sample is an integer multiple of o% — o\. The number of 
"quanta" M depends on the value of V. A typical shape of the peak for M max = 2 is shown 
in Fig. 4 (d). Each step in a^P results from the switch of some saddle-point from a "black" 
to a "white" bridge. Note that the resistance (0) does not demonstrate regular scaling 
with the aspect ratio of the sample since M is just a random integer varying from sample 
to sample. The probability distribution for different M does, however, depend on the shape 
of the sample. For instance, if r2 — r± ~ r±, the most probable values are M = 0, 1, or 2. In 
a very narrow ring, W <C r, the most probable values are close to M = C/A Q 3> 1, where C 
is some numerical factor. In the latter case, the resistance will seem to scale properly with 
the sample dimensions, as if the sample was homogeneous, except for the wrong numerical 
factor. The quantization of i?tot will be difficult to see. 

Our prediction of the longitudinal resistance quantization allows a comparison of our 
theory with the experiments of Ref. pT| . Although neither of the two samples shown in 
Fig. 3 reveals the saturation of the peak heights which we expect to occur at low enough 
temperature, still, we can see that the heights of all peaks become close to each other 
at the lowest temperature studied, T = 14mK. Hence we assume that saturation occurs 
not very far below this temperature and use the values of cr^ p obtained at 14mK as good 
approximations to the zero-temperature values. Taking account of the aspect ratios quoted 
in the caption to Fig. 3, for both samples we obtain R tot = A /a^P ~ h/e 2 within an 
accuracy of 10% (we averaged cr^ p (T = 0) over 5 values for different peaks in sample B). 
This corresponds to M = 1 in Eq. flT3"| ) which is consistent with our prediction. Thus, at 
least for these data, the "fluctuation" of a^P between samples is simply correlated with the 
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different aspect ratios, A , which were used to calculate a x ^P = A /R tot . 

IV. BOUNDARY STRIP FORMALISM 

So far, we have considered the simplest geometry of the Corbino disc. Although this 
method can directly produce the value of a xx for a homogeneous system, it also has some 
obvious disadvantages. Firstly, it does not allow one to measure the Hall conductivity. Sec- 
ondly, it leaves no hope of separating <t£* p into the contributions arising from the bulk and 
those arising from the edge. Both disadvantages stem from the fact that only one indepen- 
dent experimental parameter (the two-terminal resistance) is obtained in this method. To 
permit a larger number of independent measurements one must consider another geometry, 
such as the Hall bar. As we will show in this section, macroscopic inhomogeneities in a 
finite-sized Hall bar lead to a similar decrease of the observed peak heights in as for the 
observed peak heights of cx|* p in the Corbino disc geometry. In Sec. V, we will suggest a 
method by which this edge effect may be compensated in the Hall bar geometry, at least, 
over the temperature range for which this contribution is small. At the lowest temperatures, 
T <C T s2 , as clearly follows from analysis in Sec. Ill, all information on the bulk properties 
is lost beyond recovery. 

For the remainder of this paper we will therefore focus on the regime of small edge 
corrections. In this limit, it is possible to develop a theory in which the edge effects can be 
accounted for by a "boundary strip" . Within this formalism, the Corbino disc and Hall bar 
are thought of as homogeneous samples with an infinitesimally thin layer attached at their 
boundaries to account for the edge-effects. This "boundary strip" is characterized by an 
impedance matrix which linearly relates the currents in the strip to the potential gradients 
at the boundary. Effectively the presence of this strip changes the boundary conditions on 
the sample. Such an approach is valid provided (i) the region along the edge responsible 
for the edge effects is narrow and (ii) the edges are homogeneous along their length. The 
first condition is necessary since otherwise one has to take into account, in addition to the 
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gradients, second and higher derivatives of the electric potential to describe the boundary 
impedance. 

In physical devices, there are two main types of edge effects that are important for 
transport. Firstly, the finite-size effects due to the macroscopic inhomogeneities of the 
interior of the sample, as discussed above. This type of inhomogeneity, if present, is equally 
important for the Hall bar and the Corbino disc. Secondly, smooth-edge effects caused by a 
gradual change in the electron density when approaching the sample edge. One should expect 
the latter effect to be much stronger in the Hall bar geometry in which electron density at the 
edge is zero. In principle, the boundary strip formalism is rather general and can be used to 
account for both types of inhomogeneities. For the macroscopic inhomogeneities, both of the 
conditions outlined in the previous paragraph are satisfied if the correlation radius £t is much 
less then the sample width and, in the Hall bar, the distance between the voltage probes. 
Both conditions would also be satisfied for a Hall bar with a smooth edge if the voltage 
probes were small enough as not to affect the edge properties. In real devices, however, the 
probes are macroscopically large and interrupt the edge strip. Within the model of non-local 
resistance proposed by McEuen et al, based on the edge-state formalism, the contacts have 
a strong effect on the current distribution, since they force equilibration between different 



edge channels. Such a model is a clear example of a case in which the edge is not 
homogeneous along its length, and for which the boundary strip formalism does not apply. 
We will therefore focus only on cases for which the edge of the Hall bar represents a sharp 
cut in the inhomogeneous sample, i.e. the characteristic width of the edge region is less than 
the magnetic length. The boundary properties are then independent of whether the edge is 
to vacuum (Hall bar) or to metal (Corbino disc). Below we will show how the conductance 
properties of such an edge can be related to the parameters of the two-phase model. First, 
we will derive the properties of the boundary strip for a periodic array of the two regions, 
and then we will discuss how these are modified for a random distribution. 

As we have seen in Sec. Ill, the corners formed at each edge by alternating phases 
create an effective strip with properties distinct from that of the interior of the sample. 



While discussing the Corbino disc, the boundary strip has been characterized by a single 
parameter, the contact resistance. This is defined as the ratio of the voltage drop across 
the strip V to the current j±L crossing the strip, under the condition that the electric 
field component along the edge E\\ is zero. In general, quite different boundary conditions 
may be applied. For instance, in the standard Hall bar measurement, j± = 0, E\\ ^ 0. We 
need to develop a quantitative description for the boundary strip which does not depend on 
specific boundary conditions and which therefore applies to all geometries. We consider four 
variables: V,j±,E\\, and the additional current / flowing along the edge in the boundary 
strip. The four parameters are related by a contact impedance matrix E, as given by 
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which is analogous to the conductivity tensor in the bulk. Due to the low symmetry of 
the boundary strip, all four components of E are a priori independent. The whole in- 
homogeneous sample can be thought of as consisting of a homogeneous interior, with the 
conductivity tensor of the infinite system a* and with the same dimensions as the original 
sample, and an infinitesimally thin boundary strip described by the impedance matrix E. 

We will now determine the boundary strip impedance matrix for a simplified model: a 
periodic two-phase system in which "black" and "white" are regularly distributed as in a 
chessboard, and all vertices have identical scattering properties. Since, for a thin boundary 
strip, the sample shape is not important, we choose to study a convenient geometry, in 
which the sample is a long rectangle with rows of vertices aligned parallel to its sides. Let 
the width of the sample be W = Nd, where d is the lattice constant, and N is an integer. 
A section of a long sample is shown in Fig. 5(a). Although we are ultimately interested in 
the case of large N, the periodicity of the problem enables us to use a system with a small 
number of vertex rows (in Fig. 5(a), N = 3). Our objective is to replace this system by 
an equivalent homogeneous sample with thin boundary strips attached, Fig. 5(b). The new 
system is characterized by the electric field E°° and the current density j°° in the interior, 
and the four variables V, E\\, I,j± describing the boundary strip. The positive directions for 
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the potential drop V and for the other three variables are shown by arrows in Fig. 5(b). All 
these currents and fields are assumed to be uniform, representing one particular case which 
is sufficient to evaluate matrix (|b3|). We will consider the field E°° and the current density 
j°° to be given as boundary conditions. Even though these two vectors are related by the 
conductivity tensor of the infinite system a*, in what follows, we are not going to use this 
relation, and will treat E°° and j°° as independent variables. 

The criteria for the equivalence of the two systems are as follows, (i) The field .Eoo and 
the current density in the new sample must be the same as the average electric field and 
the average current density which would be in the original sample if the latter was infinite; 
that is, the properties of the interior are independent of the presence of the boundaries, 
(ii) The standard continuity conditions j± = j™, E\\ = E^ must be satisfied, (iii) The total 
current along the new sample 

I tot = 2I + Wj™ (15) 

and the total potential drop across the new sample 

V tot = 2V + WE™ (16) 

must be the same as the corresponding quantities in the original system. 

To apply these rules, consider the current-field distribution in the original system, Fig. 
5(a). As shown in Refs. @,0], each 4- vertex is characterized by the white-to- white current 
Ii and the black-to-black current Jj. Each arrow in Fig. 5(a) denotes a set of current lines. 
When passing through the white or black square, the current spread all over the square, then 
focuses at a corner where it can pass to another square. We remind the reader that current 
lines cannot cross the boundaries between quantized regions so that the specific distribution 
of these lines inside any square is irrelevant. Each side of a square is at constant potential. 
The potential drop between the two sides forming a (white or black) corner are given by 

= Vi=± (17) 
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where (J») is the total current focusing at the corner. Since the system is periodic, and 
average fields and currents are taken to be uniform, the set of local currents Jj, Jj must 
also be periodic. The pair of currents at a vertex takes either of two values, 1%, J\, for odd 
vertices, and I2, J2, for even vertices. The current crosses a sample boundary by focusing at 
3- vertices as shown in Fig. 5(a) for the case 02 > <J\. This occurs at every other 3- vertex 
since current lines can focus only if a xy experiences a step-like increase to the right when 
looking from the edge into the sample (see Appendix) . The splitting of arrows schematically 
shows the splitting of the sets of current lines. A diagram similar to Fig. 5(a) can be drawn 
for the equipotential lines, except the notations Ii and Jj must everywhere be replaced by 
Ui and Vi respectively. 

The components of the average current density in an infinite sample can be obtained as 
an average of the corresponding local currents over vertices of the two types, as given by 

ff = /l " /2 2 + /' + J2 - (18) 
)7 = h + h-j 1 + h (19) 

Analogously, for the electric field components we have 

where we used Eqs. (|17|) . The total current along the sample I to t can be evaluated by adding 
the local currents crossing the vertical dashed line shown in Fig. 5(a). The result can be 
written as 

I iat = NdjZ> + 2I, (22) 
I = \(h+h + Ji-J2), (23) 
where is given by Eq. fll8|) . The total potential drop across the sample works out to be 

V tot = NdE™ + 2V, (24) 

V lf!l±I> + ±Z*\ (25) 
4 V o 2 <J\ J 
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where is given by Eq. (pip . As we can see from condition (iii) [Eqs. ( |I5|) , ([IB"!) ] the 
parameters J and V introduced in Eqs. (|23| ) and ( |25|) represent, by definition, the effective 
current in the strip and the effective voltage drop at the strip, respectively. Notice that the 
right-hand sides of Eqs. (p3|), (|25| ) depend only on I\ + I 2 and J\ — J 2 . Using the system 
of Eqs. (0), (f^]), we can express these two linear combinations in terms of components 
j™, E™. On the other hand, from condition (ii) we have = j± and E™ = E\\. As a result, 
we arrive at two equations relating I and V to j± and En . This can be written in the form 



([UP with matrix £ given by 



( d . 1 \ 

--(t7 2 -«7i) -(<7 2 + <Ti) 



y -2^2 + ^1) ^(^2 - o-i) y 



(26) 



Note that during this derivation we did not make any particular assumptions about the 
conductivity tensor in the bulk a* since we did not specify components of j°° and E°° . 

Let us discuss possible modifications of our result for the case in which the two-phase 
system is random with a correlation radius £t, Eq. (0). Obviously, the Hall components 
S12 and S 2 i in Eq. (p6|) will not change since they simply reflect the fact that black and 
white phases, in the vicinity of the inter-plateau crossover, share the edge equally. Then, the 
diagonal components S 22 have to be proportional to o~ 2 — 0\ since they originate from 
the dissipative resistance of 3- vertices, Eq. (0). The lattice period d in the two components 
En, S 22 has to be replaced by some average lengths, d\ and d 2 , respectively, both of the 
order of £y. One might expect a priori that d\ and c? 2 could differ by a numerical factor, 
since directions along and across the edge are not equivalent. However, as we show below, 
both lengths are exactly equal, d\ = d 2 = d. 

We will employ a duality that exists between the current and field distributions in 2D 
conductors. |22| Let us imagine that the system in Fig. 5(a) (with many rows of vertices) 
is randomized, i.e., squares are distorted and vertices are not identical. The currents Ii, Ji 
and the voltage drops U^Vi at corners are no longer periodic. Although the system is 
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fully characterized by the discrete set of currents and voltages at the vertices, it will be 
convenient for now to consider the local current density distribution j(r) and the local 
electric field E(r). Both functions satisfy the continuity conditions that the number of 
current or potential lines entering and leaving a given (black or white) square are equal. 
Suppose that we have found the contact impedance matrix in this system which, as explained 
above, has a form 



—4(02 -(Ti) 
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1 



sPi + cr lj 



1 

do 



lcr 2 - cri 



(27) 
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Let us now map our system onto a new (primed) system with the same geometry of the 
phase distribution and with the new current density and electric field 



E\r) = [zx j(r)], f(r) = [z x E(r) 



(28) 



The steady-state conditions, V ■ j' — 0, V x E' — 0, are obviously satisfied in the new 
system, if they are satisfied in the original system, V • j — 0, V x E — 0. As follows from 
Eqs. (|28f) , the new quantized Hall conductivities in the black and white are 



-pk 



■l/a k , k = 1,2. 



(29) 



The average components I ,Vpriortopublication, j±, Eu characterizing the boundary strip 
are transformed in the same way as the components of the local current and electric field in 
Eqs. (H), that is 



(30) 



V' = I, E\ = -ju 
I' = -V, f ± = E {l . 

Since the geometry of the phase distribution in the primed system does not differ from the 
original one, the variables should be related to EL V via the same impedance matrix 
E given in Eq. Q27D, except a k should now be replaced by a' k . Using Eqs. (PD|), we 

arrive at the relation 
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(31) 



Comparing Eqs. ( Ill ) and (51), we see that P = As one can easily check, this is only 
consistent with Eq. fl27|) if c?i = d%, which proves our assertion that d\ — d^ = d even for a 
random system. 

Thus the presence of the edge is equivalent to a fictitious homogeneous anisotropic strip 
of width d/2 ~ £t with the local resistivity tensor components 



Pyx 



Pxx 



Pxy = (Pi + P2)/2, 



-P 



xy 
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yy 



(Pi-p 2 )/2. 



(32) 
(33) 



The unusual fact that the dissipative resistivity in the direction perpendicular to the strip 
is negative deserves comment. A real physical strip (or a layer if in 3D) with well-defined 
geometric boundaries cannot have a negative net diagonal resistivity in any direction, since 
this would contradict the Second Law of thermodynamics. However, this is not the case 
here: the effective contact strip in our discussion has no real geometric boundary which 
could be drawn, for instance, inside of the sample in Fig. 4(a). The strip describes small 
corrections to the net conducting properties of the sample which, as a whole, has a positive 
dissipation. 

We can now reobtain the contact resistance in the Corbino geometry, say, that from the 
inner contact ARi. Putting E\\ = 0, from Eqs. (0), (^) we have 

V d 



(34) 



j_i_2nri 27rri(0- 2 - a x ) 
which coincides with the corresponding term in -R t ot, Eq. (|TT|) , if one puts d = ^t- 

Consider now the standard Hall bar measurement in which the current flows parallel to 
the edges, j± = j{ 



jOO 

'y 



0. Relation (|31| ) taken with d\ = d 2 = d yields 
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£|l = ^p/, V = p b y J, (35) 



where p b xx and p b are given by Eqs. fl33|) and (|32|). The electric field in the interior is 
homogeneous and given by 

where p* = (a*) 1 is the resistivity tensor of an infinite sample. Expressions for the exper- 



imentally measured components of the resistivity tensor p e x x x p = WE\\j I tot , p^ p = V to t/I\ 



tut 



can easily be obtained from Eqs. fli"5D, (|TE|), fl3"5|), and The result has a form 



1 hot 1 d ,„„. 

+ 77nr> ( 37 ) 



exp * _ rxx 

"yx ~ j ~ Pyx I.; b 
tot ^P 2 a 



(p£* - • (38) 



The last expression was expanded in terms of the small parameter d/W. 

Thus, the presence of the edge leads to a negative correction in the measured p^P, as it 
did for a x x x p in the Corbino geometry. Moreover, if p 1 — p 2 pi, the relative magnitudes 
of both corrections are the same provided the Corbino disc is narrow and has the same 
ratio d/W as the Hall bar. As we have seen above, this fact is related to the current-field 
duality. Physically, the negative correction to a^P in the Corbino measurement results 
from an additional voltage drop at the edge, and the negative correction to p e xx v in a Hall 
bar measurement results from an additional current trapped at the edge. Our conclusion 
agrees with the experimental observation that the peak heights tend to decrease at low 
temperatures in both geometries. []lT|L| The correction to the measured Hall conductivity 



changes sign in the middle of the crossover region where p yx = (pi + p%)/2. As a result, 
the characteristic width of the transition for both functions p^ij 1 ) and p^(v) is somewhat 
increased by a relative factor of d/W. 

When calculating the impedance matrix above, we considered the particular case in 
which the average current density and electric field in the interior of the equivalent sample 
shown in Fig. 5(b) are homogeneous. In the following section we will discuss "non-local" 
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resistance measurements for which this is not the case. We now discuss the applicability of 
the boundary strip approach for an inhomogeneous electric field. Suppose that the electric 
field in the interior of the equivalent sample E°°(r) varies with some characteristic length 
Ie-i where lg ^> d. Then the four parameters of the boundary strip entering matrix relation 
(|H|) will also depend on the coordinate along the edge x. Let us express the current density 
in the homogeneous interior in terms of the pseudoscalar ip(r) as given by 

j(r) = [zx V^(r)], (39) 

which is always possible since V • j — 0. Analogously, the electric field can be written in 
terms of the electric potential as E°°(r) = — V<p(r). At the boundary strip, ip experiences 
a step from ip\ to ip2- The current inside the strip / and the current crossing the strip j± 
are given by 

t i i ■ I (dip! dip 2 \ , An s 
Analogously, V and En are given by 

where <pi and 02 are the potentials at the outer and inner sides of the strip, respectively. 
Thus, the current across the strip is defined as the average of the currents crossing the 
inner and the outer sides of the strip. Since the current inside the strip / depends, in the 
general case, on the coordinate x, the two currents may be different. Similarly, En is the 
average of the parallel components of the electric field at the two sides of the strip. Such 
a choice of the definitions of En and jj_ ensures that all parameters entering the relation 
(|I4|) are expressed via first derivatives of functions ip or <p\ I and V are discrete derivatives 
(differences) in y, and j± and En are continuous derivatives in x. As a result, the matrix 
relation ( p~4"D represents the correct description of the conducting properties of the edge to 
first order in the small parameter d/lE- To increase the accuracy to second order, one would 
have to write a matrix relation which also includes second derivatives of ip and <fi such as 
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dipi/dx — dip2 / dx , d 2 ip\ / dx 2 + d 2 ip2/dx 2 etc. In this and the following sections, we restrict 
ourselves to first-order effects in djl^. We note that the aforementioned equivalence between 
the effective boundary strip and a homogeneous strip with a resistivity tensor p b is also only 
correct to first-order in this parameter. 

V. COMPENSATING THE FINITE SIZE EFFECT: "NON-LOCAL" 

RESISTANCE 

The objective of a standard QHE transport experiment performed on a large sample is 
to extract the bulk conductivity tensor characterizing an infinite system. As we have argued 
above, at low temperatures comparable to T S 2, finite-size effects become noticeable. It would 
be very useful to devise a method by which these edge contributions could be separated from 
the measured resistances. In this section we show how this may be achieved in systems for 
which the boundary-strip formalism of section IV applies: that is, for samples with a sharp 
edge and for which the edge effects are not too large. 

As shown in the previous section, at temperatures which are not too low, the edge effect 
can be described by a matrix X (or P = which, for the two-phase model with a 

sharp edge, contains a single unknown parameter - the average length d ~ In order to 
determine this parameter, one additional measurement beyond the standard measurements 
of R xx and R xy in the Hall bar geometry is required. In what follows, we suggest a way 
in which d may be extracted from a measurement of the enhanced "non-local" resistance. 

Unlike the standard Hall bar measurement, for which the current passes along the sample, 
in a "non-local" measurement the current is forced to cross the sample between probes 1 and 
2 on opposite long sides of the sample, Fig. 5(b). The "non-local" resistance is determined 
from the potential difference between a second pair of probes 3 and 4, -R n ioc = V34//12. In 
a homogeneous sample, as follows from standard electrostatic considerations, -R n i oc should 
decay with the distance between the current and voltage probes, L, as given by the series 
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tfmoc = p*„ (C ie -^ W + C 2 e'^ L ' w + ...) , (42) 

where G\ ~ C2 ~ 1 are numerical coefficients determined by the shape of the contacts. At 
L > W, the resistance is dominated by the first exponential in this series. As found in Ref. 
|2l]| , the experimental value of -R n ioc observed at low temperatures is much larger than that 
predicted by Eq. (£2). This effect clearly indicates the existence of currents localized near 
the edge, in addition to the current passing in the interior of the sample. 

The most common explanation for the enhanced "non-local" effect invokes the presence 
of a smooth edge to the sample, at which the electron density vanishes slowly. This would 
lead to the appearance of one or more quantized Hall strip (s) at the edge, within which the 
scattering can be very small, even when the bulk of the sample is in the region of the peak in 
p xx . The low-dissipative strip can trap a noticeable portion of the current in the sample. 
In the edge-state transport language, this can be formulated as a poor equilibration between 
different edge channels: those at the edge of the sample and those in the bulk. |13|,^4| ^ 
is also possible for an enhanced "non-local" effect to arise, even when the edge is abrupt, 
as a result of the random macroscopic inhomogeneities which we have discussed in previous 
sections. As we have seen above, the effective boundary strip traps an additional current / 
along the edge which causes a decrease in the observed R xx . As we will show in this section, 
the same edge current causes an enhanced "non-local" effect which can be thought of as an 
effective increase in the sample width. 

We present a simple quantitative theory which allows one to relate the corrections in 
R xx to the enhancement in R n \ oc . Since our derivation is based on the phenomenological 
boundary strip description, this approach is rather general and can be used for a class of the 
edge models. Moreover, although the edge properties are described by the four components 
of the matrix p b which are expected to depend on a specific edge model, the relation between 
the corrections in R xx and the enhancement in i? n i oc turns out to be universal in the sense 
that it does not include any of these components. The edge corrections to the measured 
Hall resistance R yx can be found from the non-local effect in the same manner. In this case, 
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however, one has to know the component p b of the boundary strip matrix. 

We begin by deriving the current and potential distributions for a non-local measurement 
on the effective sample shown in Fig. 5(b). To do so, it is convenient to express the current 
density in terms of a pseudoscalar ip(r), as defined by Eq. (|39|). Correspondingly, the 
currents inside and across the strip are given by Eq. (H0|). The distribution ip(r) in the 
interior of the sample satisfies the Laplace equation 



V 2 ip = 0, 



(43) 



which follows from the conditions V x E = 0,p* x 7^ 0. In addition, the parallel component 
of the electric field in the interior E x and the electric field in the strip En must match, as 
given by 

E x {y = -W/2) = E h 

Ex = Pxx3x — PyxJyi 



E ll= p b xx 2I/d-p b yx j ± , 



(44) 



where the matrix p b is defined in Eqs. 
can be written as 



. Using Eqs. fl3"9]), (HO), the last condition 



2 b , i i s Pyx dip 
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\ Pyx 2 



dx 



+Pa 



dy 



0. 



(45) 



where all partial derivatives are evaluated at the lower edge in Fig. 5(b), y = —W/2. The 
matching condition for the upper edge is analogous. 

The solution of Eqs. (f43l), (|45f) depends on V ; i( a; ) , i- e -; on the conditions on the outer 
boundary of the sample which are set in the experiment. Let us assume that, in the "non- 
local" resistance measurements, the current 1 12 enters and leaves the sample at small contacts 
1 and 2 which are positioned at x — 0, Fig. 5(a). Then we have ipi — In/2 at x > and 
■01 = — -Z12/2 at x < 0, with a step at x = 0. We will restrict ourselves to finding ip(r) far 
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from the current probes, x 3> W. In the limit x — > oo, all currents vanish so that ip2 — ► ipi- 
We will look for an asymptotic solution of Eqs. (|43|), (|45|) in the form 



i)(r) = J 12 /2 + C cos(ky)e~ kx (46) 

which is symmetric about y = and satisfies Eq. ((4~3|). Substituting this anzatz and 
dipi/dx = into (|45"D, we obtain the following equation for the decrement k 



XX 



We have to choose the smallest solution of this last equation which tends to to n/W as 
d/W — > 0. To first order in d/W, we obtain 

k = ^~, W eS = W + d^. (48) 
W eS p xx 

Clearly, the voltage between probes 3, 4 decays the distance L in the same exponential 
way as the current density, such that -R n i oc oc exp(— 7rL/W e g). Thus, as far as non-local 
resistance measurements are concerned, macroscopic inhomogeneities make the sample ef- 
fectively wider. In the region of the peak in p* xx (v) where p* xx ~ p b xx , the effective width 
increases by of order the correlation radius £t- 



Notice that the formula (|37D for the measured diagonal resistivity can now be rewritten 



as 

Pxx Pxx ' 

(49) 

which means that effect of inhomogeneities on the "local" resistance can be thought of as 
the same increase in the sample width: R xx = p e xx p L/W = p xx L/W e g. This immediately 
gives a simple method by which the edge contribution to p exp can be eliminated to first order 
in d/W . First one must measure the "non-local" resistance -R n ioc with two pairs of voltage 
probes situated at different distances L±,L2 (larger than W) from the current probes. The 
effective width can be calculated as 



W eS = n(L 2 -L 1 ) In" 1 



-Rnioc (-^l) 
-Rnioc (-^2) 



(50) 
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Combined with a knowledge of the geometric sample width W this can then be substituted 
in Eq. ( fj9|) to obtain p* x . (Since precise measurements of the geometric width of the sample 
W may present some difficulties in experiment one may determine it in the same way as W e s, 
but at B = when weak inhomogeneities are not important.) Using the same method, we 
can extract the bulk value of the Hall resistivity p* x from observed Hall resistance R yx = p^, p . 
From Eqs. (||), fl3|) we get 



Pyx — P C yx + (Pyx P — Pyx) (^y ^ ) " ' ' ' 



It is worth emphasizing that when deriving expressions fl49| ) - (|5l|) , we did not use any 
particular values for the components of the boundary strip matrix p b . This means that 
the method of compensating the edge effect we just have suggested is rather general and is 
not restricted to the "black-and-white" regime, T -C T sl , studied in Sees. Ill and IV. For 
instance, the method can be applied just as well at higher temperatures, T > T s \, when 
a significant part of the current is carried by the intermediate "grey" region in Fig. 1(c). 
This seems to be the case in experimental data shown in Fig. 3 at high temperatures when 
a^P decreases with increasing T. Relation (^) used to compensate the edge effect in the 
diagonal resistivity is quite universal in the sense that it does not depend on any components 
of the boundary strip matrix. In order to exclude the edge effect from the Hall resistivity 
using Eq. ([51]), one has to know the specific value of p b yx . We note, however, that the 
value Py X = (pi + P2)/2 is more general than its derivation based on the "black-and-white" 
model which we gave in Sec. IV. It merely reflects the symmetry of the local Hall resistivity 
distribution p yx {f) with respect to (pi + p^jl which is conserved as long as the sample is 
effectively very inhomogeneous, i.e. 5ut ^ <5z/ . 

We emphasize again that this method applies only to first order in d/W. It yields a 
partial compensation of the edge effect at T > T s2 , and is useless at T T s2 when all 
measured parameters have already saturated in temperature. 

Our conclusions are easy to test experimentally. If the underlying model is appropriate 
in a particular sample, the application of this method should produce a wide range in the 
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temperature dependence of cr^f x (T) within which the diagrams a xx vs. a xy are close to 
the "universal semicircle". In other words, by compensating the edge effects, the maxima 
in temperature like those shown in Fig. 3 can be broadened and brought much closer 
to 0.5e 2 /h. At the lowest temperatures, when the finite-size effects dominate the transport 
properties (the "saturation" regime), the method eventually breaks down and no universality 
can be retrieved. In fact, such a behavior close to the universal prediction was already 
observed, without any special methods, in the integer regime in the temperature ranges 
T = (0.5 -2)K [Ref. 0(a)], and T = (2-4)K [Ref. [g^]. In these low-mobility samples, the 
high-temperature decrease of the peak heights starts late when the size effects are already 
small. 



VI. CONCLUSIONS 

We have shown that low-temperature measurements of the quantum Hall effect in inter- 
plateau regions are very sensitive to even weak macroscopic density inhomogeneities in the 
sample. The inhomogeneities may result in strong finite-size effects even in samples which, 
from a conventional point of view, are very large. This may be a reason why the predicted 
universal behavior of the transition regions in an infinite macroscopically homogeneous sys- 
tem is so hard to observe experimentally. Within our model of macroscopic inhomogeneities, 
we were able to account for both the decrease of the conductivity peak heights at low tem- 
peratures and the curious "non-universal" scaling of the peak heights revealed in some 
samples. In the low-temperature limit, the experimentally observed peaks in a xx in the 
Corbino geometry are shown to saturate at values proportional to the differences between 
adjacent plateaus in a xy . (Analogously, the peaks in p xx observed in the Hall bar should 
be proportional to differences in p xy .) The proportionality factor depends on the specific 
realization of disorder and fluctuates strongly between samples. The model also predicts 
unusual peak shapes which show quantized plateaus in the longitudinal conductivity. The 
experimental peak heights obtained in Ref. []TTJ at the lowest temperatures are consistent 
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with this quantization. 

Finally, we also showed that for a Hall bar with a sharp edge, there exist simple relations 
between the enhanced non-local resistance and the size corrections in R xx and R xy when 
these corrections are small. These relations can be used to separate the edge effects from 
the bulk tensor components p xx and p xy . 
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APPENDIX: 

EFFECTIVE RESISTANCE OF A FORK- VERTEX 

In this Appendix, we calculate the effective resistance R 3 of a fork vertex, a relatively 
simple example of which is shown in Fig. 2b. The current distribution in the vertex is shown 
by the arrows which denote beams of current lines. The beams split and focus at "simple" 3- 
and 4- vertices which form the fork vertex. (Recall that, since the longitudinal conductivity 
is assumed to be vanishingly small in black and white regions, the current lines cannot cross 
the border between phases or the metal boundary other then at these simple vertices.) The 
effective resistance R 3 can be found from the net Joule heat Q = R3I 2 , where I is the net 
current entering the metal contact. Although the problem would appear quite complex, it 
turns out that the dissipated heat does not depend on the specific structure of the vertex. 
We will now show this from very general arguments. 
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Let us draw an imaginary circle (the dashed line in Fig. 2b) enclosing all of the structure 
of the vertex. Consider the total currents 7 1; 7 2? and I crossing this circle in the black, white, 
and metallic areas, respectively. From the current continuity condition, I = I2 — I\, and 
from the condition that electric potentials at points 1, 3 of the metal boundary must be 
equal, I\jo\ = I2/&2, we find 

/i = -^-J, I* = -^-I. (Al) 

cr 2 — <j\ 02 — 0\ 

that is, only one of the three currents is an independent quantity. 

The total Joule heat inside of the circle is determined by the general expression 

Q = J E(r)j(r)d 2 r, (A2) 

where E(r) and j(r) are the local electric field and the local current density, respectively, 
at point r inside of the circle. Since V • j = 0, V x E = 0, we can express j in terms of the 
pseudoscalar if), and E in terms of the electric potential 0, as given by 

j(r) = [zx Vif)(r)), E(r) = -V0(r). (A3) 

In this notation, the three currents /, I\, I2 can be written as 

I = if>i - ip 3 , h=1p2~ if>i, h = if)2- ^3- (A4) 

Substituting Eqs. ( |A^ ) into (|A2|) , changing order in the mixed product, and integrating by 
parts we get 

Q = - fiP{V(f>-dl), (A5) 

where the integral is taken along the closed circle in Fig. 2b. The circle can be broken into 
three segments: 3-2, 2-1, and 1-3. Within each of the segments we have V0 = (l/a^Vip, 
where the Hall conductivity ex, is equal to 02,01, and 00 (metal), respectively. Integration 
within separate segments yields 
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Using Eqs. and (|A4|) , we obtain finally 

Q = 9 , /2 x = Rsl 2 , (A7) 
2(cr 2 - (Ji) 

which yields the formula (|S|) quoted in the main text. 

Note, that if the colors in Fig. 2b are interchanged, which corresponds to interchanging 
02 and <j\ in Eq. (|A7|) , the Joule heat will be negative (recall that we assume <j 2 > <J\ 
everywhere in the paper). Since this would contradict the Second Law of thermodynamics, 
for such a vertex, all currents must be zero. Hence, only the vertices for which the rightmost 
color, as shown in Fig. 2b, is white, can be "active" (can participate in the current transfer). 

To illustrate our result (|9]), consider the simplest (without branching) example of a 3- 
vertex presented by a corner of a rectangular homogeneous conductor with the Hall conduc- 
tivity 02. Two metallic probes are attached at the bottom and at the top of the rectangle, 
the sample on the left and on the right bordering to vacuum, a\ = 0. If 2 > 0, the current 
will leave and enter metallic probes focusing at the lower left and the upper right corners 
of the sample which represent the active 3-vertices. The two-terminal resistance of such a 
sample is equal, as easy to see, to the inverse Hall conductivity I/02. This amounts to the 
effective resistance R3 = 1/202 per each active 3- vertex, in agreement with Eq. (Rp. 

As is easy to check, the same rule of selection of active vertices applies if the metal in 
Fig. 2b is replaced by a vacuum, as is appropriate for a Hall bar with an abrupt edge (Sec. 
IV). In contrast to the edge to metal discussed above, we now have zero current out of the 
edge, 1 = 0, and a non-zero voltage difference V = 03 — <j>\. The expression for the Joule 
heat Q can be obtained in almost the same manner as Eq. ([S^]), except it is now convenient 
to use, instead of Eq. (|A5|), an equivalent formula 



Q = J 0(W • dl). (A8) 

The final answer has the form 

V 2 

2(0! -0 2 ) 

which has the same sign as the right-hand side of Eq. 
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FIGURES 

FIG. 1. (a) Local conductivity tensor components a xx and a xy versus the local filling factor 
v. (b) Macroscopic fluctuations of the local filling factor, (c) Conductivity tensor distribution in 
inhomogeneous system at low temperatures. Black and white regions correspond to quantized Hall 
regions with a xy = o\ and 02, respectively, and o xx « 0. The grey color which corresponds to the 
grey strip in (b) shows the intermediate (non-quantized) region. 

FIG. 2. (a) Conductivity distribution in an inhomogeneous Corbino disc. Black and white are 
quantized regions as in Fig. 1(c), the intermediate grey region is not shown, (b) Complex "fork 
vertex" from the small square in (a) shown in magnified view. The arrows schematically show 
currents passing between different regions by focusing at simple vertices. 

FIG. 3. The temperature dependence of peak values of a xx for different transitions between 
adjacent IQHE states obtained by Rokhinson et al. Different symbols correspond to transi- 
tions: 5-6 (y), 6-7 (o), 7-8 (•), 10-11 (x), 13-14 (A). Some transitions are omitted by authors 
of Ref. |0] for sake of clarity. The geometric aspect ratios of the samples are A = 0.21 (sample A) 
and A = 0.32 (sample B). 

FIG. 4. Two basic configurations of the conductivity distribution in the Corbino disc in the 
saturation regime (T <C T S 2) shown at v\ < V < v^, where v\ and i>2 are filling factors of the 
two saddle points, (a) Configuration with a finite two-terminal resistance. Thin lines are the lines 
of equal potential, (b) Configuration with infinite two-terminal resistance. No current passing 
between contacts, (c) Predicted observed diagonal conductance a x ^ vs. average filling factor V for 
configuration (a), (d) Analogous dependence for a configuration controlled by four saddle points. 

FIG. 5. Diagram illustrating the boundary impedance matrix formalism, (a) Chessboard model 
of the two-phase distribution in a long rectangular sample. The arrows schematically show currents 
passing from one (white or black) region to another by focusing at vertices. Wide arrows show 
currents flowing from (to) metallic contacts, (b) Equivalent sample formed by infinitely thin 
boundary strips attached to a homogeneous interior of the same width as that of the sample in (a) . 
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